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High-order inertial phase shifts are calculated for time-domain atom interferometers. We obtain 
closed-form analytic expressions for these shifts in accelerometer, gyroscope, optical clock and photon 
recoil measurement configurations. Our analysis includes Coriolis, centrifugal, gravitational, and 
gravity gradient-induced forces. We identify new shifts which arise at levels relevant to current and 
planned experiments. 

PACS numbers: PACS numbers: 39.20.+Q, 32.80.Pj, 03.75.Dg, 04.80.-Y 



Atom interferometric measurements have growing ap- 
plications in basic and applied science. For example, 
atom interferometric techniques have been recently used 
to measure rotations 0, gravity gradients ^j, h/mcs U 
and accelerations Q with unprecedented precision. Fu- 
ture applications range from tests of General Relativity 
to the development of next generation inertial naviga- 
tion systems. The accuracy of atom interferometric in- 
struments, however, hinges on the accuracy of the theory 
used to connect the measured interferometric phase shift 
to the physically relevant quantities. 

In this paper, we develop analytic expressions for the 
response of commonly used atom interferometer measure- 
ment configurations to experimentally relevant combina- 
tions of rotation, gravity and gravity gradient-induced 
forces. We identify new classes of high-order phase shifts 
which are observable in current experiments, and which 
seem to be of vital importance for proposed future ex- 
periments. In particular, we show below that current 
acceleration measurements need to be corrected at the 
part per billion (ppb) level for high-order rotation ef- 
fects, that measurements of the photon recoil need a 10 
ppb correction for gravity gradients, and that next gen- 
eration optical time standards will need gravity gradient 
corrections to achieve fractional frequency accuracies of 
8v/v below 10 -17 . Related theoretical work in this area 
has treated simple (one-dimensional) models [[|, || or pro- 
vided general frameworks ||. 

For simplicity we confine our discussion to the case of 
time-domain light-pulse atom interferometers 0. How- 
ever, with minor modification our results can be applied 
to other de Broglie wave interferometry approaches [^). 
For light-pulse atom interferometers, coherent division, 
redirection and recombination of atomic wavepackets is 
accomplished via momentum exchange with an external 
driving laser field Q. Complete descriptions of experi- 
mental realizations based on this principle are given in 
Ref. In brief, an optical pulse of area 7r/2, coupling 
two stable internal atomic states |1) and |2), serves to co- 
herently divide an incident atomic wavepacket: an atom 
initially in state |1) is driven into a coherent superpo- 
sition of internal states |1) and |2), with the momenta 
of the wavepackets associated with these states differ- 
ing by the momentum of the photon used to drive the 
transition. Pulses of area it induce stimulated transi- 
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FIG. 1: Schematic illustration of the interferometer ac- 
celerometer geometry. For this realization, A<^ pl . op = 
(Sci.ACBi - S c i,ADB 2 )/fi with Sc^acbj and S^adbs represent- 
ing the classical action along the interfering paths ACBi and 
ADB 2 . 

tions which exchange the internal states associated with 
the atomic wavepackets, while changing their momenta 
by the photon recoil momentum. Atom interferometers 
are realized from sequences of these pulses. For exam- 
ple, 7r/2 — 7r — 7r/2 pulse sequences have been used to 
build atomic interferometers analogous to optical Mach- 
Zehnder interferometers, as illustrated in Fig. [I] |t]]. 

The total phase shift between interfering paths, 
A^totai, can be broken into three contributions: 

A</>total = A0 prop + A0i ascr + A0 sop , 

where A(f> prop is the phase shift due to wavepacket prop- 
agation between the interrogating optical pulses, A</>i ascr 
is the shift acquired during the laser-atom interactions 
used to manipulate the atomic wavepackets, and A0 sop 
is the shift due to the (possible) final spatial separation 
of the interfering wavepackets at the interferometer out- 
put port |lC| . Our approach neglects terms originating 
in the spatial extent of the wavepackets || and invokes 
the short-pulse (high Rabi frequency) limit for the opti- 
cal interactions in order to clarify contributions arising 
from inertial forces. We briefly summarize expressions 
for A^ prop , A0i asor and A4> scp below. 
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The A0p rop term is obtained using the Feynman path 
integral approach jy], which involves calculation of the 
difference between the classical actions associated with 
the interfering wavepacket trajectories. This is illus- 
trated in Fig. 1 for the case of the accelerometer im- 
plementation. The actions S are obtained by integrating 
the Lagrangian L over the classical trajectories T asso- 
ciated with the mean positions of each wavepacket, e.g. 
S = / r L(r(t),v(t))dt, with the path T determined by the 
classical trajectory (r(i), v(t)). This approach is formally 
correct when L is at most second order in position r(t) 
and velocity v(f) (see Ref. We find the classical 

trajectories T through integration of the Euler-Lagrange 
equations for L. 

The phase difference resulting from laser interactions 
A'/'iaser is discussed in detail in Ref. This shift is 

obtained from solution of the Schrodinger equation for a 
resonantly driven two-level quantum system. In the limit 
of impulse excitations, this solution reduces to the follow- 
ing rules for the evolution of the probability amplitudes 
of states |1) and \2): 



|l)-»iexp[ik-r(<p)]|2); \2) 



«exp[— ik • r(£ p )]|l). 

(1) 



Here k is the propagation vector associated with the cou- 
pling laser field(s) and r(t p ) is the mean position of the 
atomic wavepacket at the interaction time t p . Time de- 
pendent terms associated with the frequency of the driv- 
ing fields have been suppressed (as they cancel for the 
time symmetric pulse sequences considered here), and we 
assume the initial laser phase to be constant during the 
interferometer sequence. As an example, for the pulse 
sequence of Fig. [j], the resulting laser-induced phase dif- 
ference for an ensemble initially prepared in state |1) and 
detected in state |1) is: 



A(j> 



laser 



k • (r A - r c - r D + r B J. 



(2) 



The shift Acf> scp arises when the classical positions of 
the two interfering trajectories do not coincide at the 
exit beamsplitter (as illustrated in Fig. |l|). This shift is 
adequately approximated by 



A</> scp = p • Ar/h, 



(3) 



where p is the mean momentum of the wavepackets in 
a given output port, and Ar is the spatial separation 
between the centers of each wavepacket at the time of 
the last optical pulse |L3]]. In certain special cases - for 
example, that of uniform gravity - this contribution is 
zero. 

The above phase shifts were evaluated in a geocentric 
reference frame fixed to the surface of the Earth, as shown 
in Fig. |^. In this frame, the Lagrangian can be written 
as: 



L(r,v) 



T +g-r 



-nTijTj + n- ((r + R) x v) 



+ i(fix (r + R)) 2 



(4) 
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FIG. 2: Reference frame for the calculations. The z-axis is 
chosen to point away from the Earth center and the x-axis is 
tangent to a latitude circle pointing east. 



where to is the atomic mass, R is the displacement from 
the center of the Earth, f2 the Earth rotation rate, g the 
acceleration due to gravity and Xy is the gravity gradient 
tensor. 

Analytic solutions were obtained in a perturbative ap- 
proach [jlO|, where the action was evaluated for the full 
Lagrangian L over trajectories determined by the La- 
grangian L which neglected gravity gradients, it e.g.: 



L(r, v) = to y— + g • r + fl ■ ((r + R) x v) 
+ 1 - (fix (r + R)) 2 



(5) 



This approximation was validated numerically as dis- 
cussed below. 

We evaluated accelerometer [|, [?J, gyroscope O, pho- 
ton recoil , and optical clock interferometer ]9| con- 
figurations, as summarized in Tables @-|v|. Each con- 
figuration is discussed briefly below. We tabulate only 
the most significant terms in Taylor expansions of the 
total phase shifts. In order to give a quantitative mea- 
sure for the different phase contributions, the tabulated 
expressions were evaluated for parameters corresponding 
to current experiments. In particular, m = 2.21 • 10 -25 kg 
(Cs atomic mass), A e ff = 426 nm (effective wavelength for 
two-photon Raman transitions), latitude 41°, Earth ra- 
dius R = 6.72 ■ 10 6 m, gravity g z = —9.8 m/sec 2 and 
gravity gradient 2T XX = 2T yy — —T zz — 2g z /R with all 
other gravity gradient terms being zero (for a spherically 
symmetric Earth). The coordinate system origin is taken 
as the mean wavepacket position at the time of the first 
optical pulse, and the phase shift is determined just after 
the final pulse. 

Gravimeter. Phase shifts for an atomic fountain gravime- 



3 



TABLE I: Gravimeter phase shifts. 



Term 
k z T 2 Q. 2 y R 

J-k T 4 n T 

\2 z yz-t-zz 

-3k z T 3 v z n 2 y 
-\k z T^g z Sl 2 

2m Kz± , lzz 

-\k z T^%R 
-^klT^ 



Phase (rad) 
-2.32 • 10 Y 
4.44 ■ 10 4 
1.08 ■ 10 1 
-6.32 
-3.11 • 10 -2 
1.81 ■ 10" 
1.21 • 10" 
9.71 • 10" 
-3.47- 10" 5 
-2.79 • 10" 5 
-2.62 • 10~ 5 



Relative phase 



TABLE II: Phase terms for an optical clock with light pulses 
parallel to gravity. 



1-2 
-2 
-i— 3 



1.9 
4.7 
2.7 
1.3 
7.8 ■ 
5.2 ■ 
4.2 • 
1.5 ■ 
1.2 ■ 
1.1 ■ 



1.0 
10 
10 
10 
10 
10" 
10" 
10" 
10" 
10" 
10" 



ter/accelerometer, based ona7r/2 — 7r — 7r/2 sequence of 
vertically propagating optical pulses, are presented in Ta- 
ble [|. For the numerical estimates we used a time T = 0.4 
s between the pulses (corresponding to a fountain height 
of 78 cm) and a vertical launch velocity v z = —g z T. Note 
that in a satellite microgravity experiment the first two 
terms would cancel, such that the gravity gradient contri- 
bution would dominate. Recent terrestrial experiments 
have demonstrated the capability to resolve phase shifts 
at or below the 10 -9 g level [[|. 

Optical clock. The optical Ramsey pulse sequence, tt/2 "f 
— 7r/2 I — 7r/2 I —tt/2 I (arrows indicate propagation di- 
rections of the interrogating pulses) is presented in Table 
Results are shown for pulses propagating along a ver- 
tical axis, with T = 0.4 s between the first and second and 
the third and fourth pulse, infinitesimal time between 
the second and third pulse and initial vertical velocity 
v z = —g z T. Atomic clocks based on this sequence offer 
the prospect of pushing the accuracy of the definition of 
the second below the 10 -18 level and are a longstand- 
ing goal in frequency metrology. For an operational time 
standard, terms linear in k z can be suppressed by pulse 
reversal techniques Q (k z — » — k z , for example), leaving 
terms quadratic in k z as possible systematic shifts. The 
largest of these tabulated is the well-known recoil-shift. 
The smaller term, which is linear in T zz , depends on the 
location of the measurement, and represents a possible 
systematic offset at the 8v/v ~ 10~ 17 level. This term 
exists for both horizontal and vertical interrogation ge- 
ometries. 







Photon recoil measurement. Chu and coworkers have 
shown that a modified form of the optical Ramsey 
method can be used for precise determination of the 
quantity h/mc s ||. The modification involves insertion 
of a series of N — 1 n pulses, of alternating propagation 
direction, between the second and third 7r/2 pulses, and 
has the effect of enhancing the photon recoil phase shift 
terms. Following Refs. || we calculate the phase 
difference between the two possible closed interferometer 
branches (see Table [II). The phase terms are evaluated 



Term 


Phase (rad) 


Relative phase 


k z T 2 g z 
k z T 2 Q. 2 y R 
-k 2 T-^- 

,? m 

k z T 3 v z T zz 
-^k, z T QzTzz 
-3k z T 3 v z D. 2 y 

h h.2rp3rp 
l^ Kz± ±ZZ 


-2.32 • 10 Y 
4.44 ■ 10 4 
-4.16 ■ 10 4 
1.08- 10 1 
-6.32 
-3.11 ■ 10~ 2 
1.81 ■ 10 -2 
1.21 • 10~ 2 
6.48 ■ 10" 3 


1.0 
1.9 ■ 10~ 3 
1.8 • 10" 3 
4.7 ■ 10 -7 

2.7 • 10" 7 
1.3 ■ 10" 9 

7.8 ■ 10" 10 
5.2 ■ 10"" 10 
2.7 ■ 10" 10 


TABLE III: Phase terms for a photon recoil measurement. 


Term 


Phase (rad) 


Relative phase 


21\ll U ir P 

m z 

Nil 7.2 rpSrp 
J^ Kz± 1ZZ 
N 2 hl2rp2rp rp 

2m h ' z± ^rec^zz 
(2JV 3 +JV)ft ^rrrrH rp 


8.39 • 10 b 
6.89 ■ 10" 3 
8.22 ■ 10~ 4 
4.36 • 10~ 5 


1.0 

8.2 ■ 10" 9 
9.8 ■ lO" 10 
5.2 ■ 10 -11 



between the first pair as well as the second pair of 7r/2 
pulses, Tree = (1/3000) s between the second 7r/2 pulse 
and the first ir pulse, T rcc = (1/3000) s between each of 
the subsequent N — 1 it pulses, and N — 1 = 30 ir pulses. 
The second correction term due to gravity gradients is 
at the anticipated level of precision achieved in recent 
h/mcs measurements. 

Gyroscope. A tt/2 — tt — 7r/2 sequence with optical prop- 
agation vectors nominally perpendicular to the mean 
atomic velocity can be viewed as a Sagnac-type rota- 
tion sensor || . In Table IV we estimate the phase shifts 
for a time-domain interferometer with parameters which 
correspond to the precision gyroscope of Ref. |Q. In 
particular, the light pulses are chosen to propagate hor- 
izontally in the west-east direction, with the atomic ve- 
locity v y = 290 m/s in the north-south direction and 
T = 1/290 s time between pulses (corresponding to a 1 m 
spatial separation of the laser interaction regions). The 
largest correction to the well known Sagnac shift (the 
leading term) is near the resolution limit of current in- 
struments, and is compensated by atomic beam reversal 
techniques. 



It is interesting to compare the above results with those 
from a perturbative treatment of gravity gradients and 



TABLE IV: Phase terms for a Sagnac rotation sensor. 



Term 



Phase (rad) 



Relative phase 



using the following parameters (chosen to correspond to 
the experiment in Ref. (l5|): T — 0.13 s for the time 



-L z^y 

2 fax 'J- 
-2k x T 3 fl y R 
-2k x T 3 Q y Q 2 z R 

h 1.2 rp3rp 



y9z 

3 ; 



4.69 
6.28 ■ 10~ 4 
-1.20 ■ 10~ 6 
-9.09 ■ 10~ 7 
3.11 • 10~ 9 



1.0 

1.3 • 10 -4 
2.6 ■ 10 -7 
1.9- 10"" 7 
6.6 ■ 10 



10 
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TABLE V: Phase terms for the accelerometer sequence de- 
rived treating rotations and gravity gradients as perturba- 
tions. 



Term 



k z T 2 g z 
k z T 2 Q, 2 y R 
k z T 3 v z T zz 
Y^k z T QzTzz 
k z T 3 v z Q. 2 y 

h t.2rp3rp 

-T7)k z T g z £l y 



Phase (rad) 



-2.32 ■ 10 
4.44 • 10 4 
1.08 • 10 1 

-6.32 
1.04 • 10 _: 
9.71 ■ 10~ : 
-6.05 ■ 10" 



Relative phase 



1.9 
4.7 
2.7 
4.5 • 
4.2 ■ 
2.6- 



1.0 
10 
10 
10 
10" 
10" 
10" 



rotations, i.e. using the Lagrangian L — m (v 2 /2 + g • r) 
to determine the classical trajectories, but evaluating the 
action with respect to the full Lagrangian L above. As an 
example, Table [v| shows these terms for the gravimeter 
configuration. Comparison with Table [i] indicates that 
the phase shift error associated with this commonly used 
approximation |l(| is ~ 5 • 10~ 3 rad or ~ 2 ■ 10~ 10 g, 
which might be resolved by current experiments. In con- 



trast, using L to estimate the classical trajectories, we 
obtain agreement at the ~ 1 /xrad level between our an- 
alytic calculations and numeric calculations which use 
the full Lagrangian L to determine the classical trajecto- 
ries, (ir&d level agreement is also obtained for the optical 
clock (Table ||) and photon recoil (Table II] ) configura- 
tions, while 200 picorad agreement is obtained for the 
gyroscope configuration (Table [V). 



In conclusion, we have analyzed accelerometer, gyro- 
scope, photon recoil and optical clock interferometer con- 
figurations. We have identified new phase shift terms 
which have their origin in cross-couplings between rota- 
tion, acceleration and gravity gradient perturbations on 
wavepacket motion. 
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